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Abstract 



We discuss Bose-Einstein correlation functions and momentum distributions 
O ■ which describe boson-emitting systems containing a central part surrounded 

i-G ■ by a large halo. If the characteristic size of the halo is sufficiently large, the 



contributions of central part and the halo to the invariant momentum distribu- 



5— i 

tion are shown to be separable. The momentum-dependence of the intercept 

parameter of the correlation function plays a central role. Almost all high 

energy reactions including lepton-lepton, lepton hadron, hadron-hadron and 

nuclear reactions are interpretable as boson emitting systems with large halo. 

The results are applied to high energy heavy ion data taken at CERN SPS. 

New insights are obtained for the description of low transverse momentum 

enhancement of pions. 
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Introduction. Hanbury-Brown - Twiss correlations were discovered 40 years ago, revealing 
information about the angular diameters of distant stars |I|. The method, also referred to 
as intensity interferometry, extracts information from the quantum statistical correlation 
function for (partially) chaotic fields. The method has extensively been applied to the study 
of the freeze-out geometry in high energy nucleus-nucleus collisions as well as in elementary 
particle reactions. For recent reviews see refs. [^],[3[ . 

Accumulating evidence indicates that the space-time structure of the pion emission in 
heavy ion reactions at the 200 AGeV bombarding energy region at CERN SPS has a peculiar 
feature, namely that the boson emission can be approximately divided into two parts: the 
centre and the halo [§-0]. The central part corresponds to a direct production mechanism e.g. 
hydrodynamic evolution or particle production from excited strings, followed by subsequent 
re-scattering of the particles. It is surrounded by a pionic halo which corresponds to pions 
emitted from the decay of long-lived hadronic resonances, like u, rj, rj' and K°, which have a 
mean decay length of more than 20 fm. These resonances give rise to mainly low transverse 
momentum pions. The kaonic halo is very rare, created by the weak decays of charmed 
mesons. The halo for 7 production is more populated, since some photons are emitted from 
the decays of abundant long lived resonances like ir°, rj and rj'. 

Within the halo of long lived resonances, the momentum distribution of the emitted 
bosons is independent from the (approximate) position of the decay. This in turn is smeared 
out over a large region due to the long mean decay length. This mechanism will result in 
bosons of similar momentum from a very large volume. This is in contrast to the core which 
when expanding fast in a hydrodynamic manner will have pions of similar momentum from 
a relatively small volume only. 

The key point is the following: Let us consider an ensemble of long-lived resonances 
with similar momentum, emitted from a given small volume of the core. The momentum 
distribution of the decay products of these resonances will be similar to each other, inde- 
pendently of the approximate position of the decay. Now the approximate position of the 
decay is randomly distributed along the line of the resonance propagation with the weight 



P(t) oc exp(— m res T res t/E res ). Thus the decay products will be emitted with the same 
momentum distribution from a volume which is elongated along the line of resonance prop- 
agation, given by Vd eC ay — Aq \ p re .s \ /(irires^res), where Aq is the initial transverse size of 
the surface through which the resonances are emitted with a momentum p res approximately 
at the time of the decay of the core, t s . 

Note that the volume of this elongated tube which emits resonance decay products with 
certain momentum is not only proportional to l/r res but it depends also on the momentum 
of the resonance in form of a factor of | p | /m. Thus the considered volume is large only 
if the decaying resonances satisfy | p | /(mr res ) >> 1 fm/c. This in turn implies that 
the simple picture that the halo can be identified with the decay products of long-lived 
resonances may need further corrections for pions at rapidity y = with low transverse 
momentum. 

Note that there is a gap in the life-time distribution of the abundant hadronic resonances 
which decay to pions: 1/T p ~ 1.3 fm/c, l/T^y* — 0.56 fm/c, 1/Fa — 1-6 fm/c and the rela- 
tively rare K* has a life-time of 1/T k* — 3.9 fm/c, while each of the long-lived resonances 
has a life-time bigger than 20 fm/c: l/r w = 23.5 fm/c, 1/1^ = 986.5 fm/c, 1/V V = 164400 
fm/c, 1/rVo = 2.76 cm/c. Thus the life-times of the long-lived hadronic resonances are at 
least a factor of 5 - 10 longer than the life-times of the short-lived hadronic resonances. The 
life-times of short-lived resonances are of the same order of magnitude than the time-scales 
for re-scattering at the time of the last hadronic interactions. Thus the decay-products of 
the short-lived resonances will mainly contribute to the core, which is resolvable by BEC 
measurements, while the decay-products of long-lived hadronic resonances will mainly be- 
long to the halo, re-defined alternatively as the part of the emission function which is not 
resolvable in a given Bose-Einstein measurement. 

In the present paper we investigate the phenomenological consequences of such a struc- 
ture where there is a characteristic length-scale of the core which is smaller than any length- 
scale that can be connected to the surrounding halo. Halo length-scales larger than 20 fm 
give rise to a sharp peak of the correlation function in the Q < %/R ~ 10 MeV/c region, 



where Q stands for the measured component of the relative momentum. This region is 
strongly influenced by the Coulomb interaction thus the accuracy of the correlation function 
determination will strongly depend on the Coulomb corrections which constitute a problem 
in itself ||. Furthermore, this Q < 10 MeV/c region is very difficult to measure due to 
the circumstance that we must properly determine two very close tracks which makes ex- 
perimental systematic errors the largest in this region of the correlation function even for 
high resolution experiments like the NA35 and NA44 high energy heavy ion experiments at 
CERN. 

We will show in this paper that even if we cannot measure reliably the correlation func- 
tion below Q < 10 MeV/c important new information about the particle production in a 
core/halo scenario can be given by analyzing both the invariant momentum distribution of 
single particles and the correlation function for Q > 10 MeV/c. We will apply our approach 
to NA44 data on S + Pb reactions at 200 AGeV/c. 

In the present paper we investigate the core/halo scenario utilizing analytical results in 
the Wigner-f unction formalism. This formalism is well suited to describe both the invariant 
momentum distribution (IMD) and the Bose-Einstein correlation function (BECF). We uti- 
lize the version of the formalism discussed in refs. [Q,[| and applied to analytic calculations 
in refs. []10|-|15|] recently. 

Review of Wigner-function formalism. In the Wigner-function formalism, the one-boson 
emission is characterized by the emission function, S(x,p), which can be considered as 



the time derivative of the (non-relativistic) single-particle Wigner function ]T(J or as the 
covariant Wigner-transform of the source density matrix, H,|3|JT6[|. Here x = (£, r) = 
(t,r x ,r y ,r z ) denotes the four-vector in space-time, and p = (E,p) = (E,p x ,p y ,p z ) stands 



for the four-momentum of the on-shell particles with mass m = \/E 2 — p 2 . The transverse 

/ (e + Pz \ 

mass is defined as m t = Jm 2 + p 2 + p 2 and y = 0.5 log stands for the rapidity. 

\E-p z/ 
For chaotic sources the two-particle Wigner-function can be expressed in terms of the 

symmetrized products of the single-particle Wigner functions. Based on this property, the 



two-particle Bose-Einstein correlation functions are determined from the single-particle emis- 
sion function solely. These relations are especially simple in terms of the Fourier-transformed 
emission function, 

S(Ak, K)= [ d 4 x S(x, K) exp(tAk ■ x), (1) 

where 

Pi +P2 

Ak= Pl -p 2 , K = (2) 

2 

and AA; • x stands for the inner product of the four-vectors. 

The inclusive /-particle invariant momentum distribution of the emitted particles reads 

as 

I da 
JVi(pi, -., Pi) = —E 1 ...E l , (3) 

a i n dp!...dpi 

where Oi n is the total inelastic cross section. These inclusive /-particle momentum distribu- 
tions are normalized to the Z-th factorial moments of the multiplicity distribution as 

(n(n-l)...(n-J + l))= f ... [*&...** Nfa, ..., Pl ) (4) 



Ex Ex 

These distributions are basically the IMD-s of particle numbers (IMD-N). 

The absolute normalization of the IMD-s (or the corresponding emission functions) is to 
be considered carefully. The IMD-s of the detection probability (IMD-P-s) 

1 En da 
Pi(Pi) = L — , (5) 

(n) a in dpi 

1 EiE 2 d6(T 
A( P i, P2 ) = — , (6) 

(n(n - 1)) a in dp x dp 2 

are also frequently used in the literature, (e.g. see refs. P, [Tl"| , |T2| | ) , which are normalized to 
unity as 



, dpi 

/— Pi(pi) =1, (7) 

E l 

dpi dp 2 

P 2 (pi,p 2 ) =1- (8) 

Ei E 2 

For cylindrically symmetric systems, the probability distribution in y and m t variables, 
d 2 n/dy/d(m 2 ) shall also be utilized. These distributions are integrals of the -Pi(p) distri- 
bution over the polar angle <fi in the transverse plane. Obviously, one has d 2 n/dy/d(mf) = 
irPi(p(y,mt,4>o)), where O is an arbitrary polar angle and Pi is independent of 4> . 

In the versions of the Wigner-function formalism presented in refs. |9,11,12], the phase- 



space integral of the emission function was normalized to 1. This is to be contrasted to the 
versions of the hydrodynamic descriptions of the IMD-s and BECF-s in refs. |T^JT5|Ji7|JT8 



where the phase-space integral of the emission function is normalized to (n). For the sake 
of clarity, we distinguish these emission functions by the subscripts w (for Wigner) and hyd 
for hydrodynamics. Obviously, their relation is given by 

Shyd(x,p) = (n)S w (x,p). (9) 

The one-particle invariant momentum distributions (IMD-P-s and IMD-N-s) may be ex- 
pressed with the emission function as 

Pi(p) = S w (&k = 0,K = p), (10) 

Ni(p) = S hyd (Ak = 0,K = p). (11) 

The two-particle BECF-s are expressed in terms of our auxiliary function, S w (Ak, K) as 

{n) 2 iV 2 (pi,p 2 ) P 2 (pi,p 2 ) 
C(K,Ak) = = (12) 

(n(n - 1)) 7Vi( Pl ) 7Vi(p 2 ) P^p^P^p,) 

\S w {Ak,K) | 2 \S w (Ak,K) | 2 
= 1 + ~1 + , (13) 

^(0,pi)^(0,p 2 ) \§ W {0,K) | 2 



as was presented e.g. in ref. |^,[|[ll|]. The corrections to the last approximation are known 



to be negligible [|llj . Effects of final state Coulomb and Yukawa interactions are neglected 
as implicitly assumed by the above equations, and completely chaotic emission is assumed. 
The formulas recollected up to this point shall be utilized in the subsequent part of this 
paper. 

A core and halo model. The Wigner-f unctions are in general complex valued functions, 
being the quantum analogue of the classical phase-space distribution functions. For the pur- 
pose of Monte-Carlo simulations of Bose-Einstein correlation functions the off-shell Wigner- 
functions are approximated by the on-shell classical emission functions obtained in a simula- 
tion P,|19|1; in case of analytic calculations they are usually modeled by some positive valued 



functions ||iO|-|T5|l . In this paper, neither of these approximations seems to be necessary. In 
fact we need only the following simplifying assumptions: 

Assumption 1. The bosons are emitted either from a central core or from the surrounding 
halo. The respective emission functions are indicated by S c (x,p) and Sh(x,p), and f c indi- 
cates the fraction of the bosons emitted from the central part. According to this assumption, 
the complete emission function can be written as 

S w (x,p) = f c S c (x,p) + (1 - f c )S h (x,p). (14) 

Both the emission function of the centre and that of the halo are normalized to unity here, 
accordingly to the Wigner-f unction normalization. 

Assumption 2. We assume that the emission function Sh{x,p), which characterizes the 
halo, changes on a scale Rh which is larger than R max ~ h/Q min , the maximum length-scale 
resolvable by the intensity interferometry microscope. In our considerations, Q m i n is taken 
as 10 MeV/c. However, the smaller central part of size R c is assumed to be resolvable, 

Rh > Rmax > R c - (15) 

This inequality is assumed to be satisfied by all characteristic scales in the halo and in the 
central part, e.g. in case the side, out or longitudinal components P0"| , |TT| of the correlation 
function are not identical. 



As mentioned in the introduction, the halo of long-lived resonances is characterized 
by large regions of homogeneity |I7|JI5|] according to this Assumption. This is motivated 



by the fact that particles emerging from decays of long-lived resonances are emitted with a 
momentum distribution independent of the approximate position of the decay, which in turn 
is smeared out over a large volume being proportional to the decay-time of the (long-lived) 
resonances. 

Note, that Assumption 2 relies on 

Assumption 0: The core does not have a self-similar, no-scale structure. The length- 
scale^) of the core is (are) resolvable by the HBT microscope. 

Based on Assumption 1 the IMD-s can be re-written as 

Pi(p) = / c Pi, c (p) + (1 - f c ) Pi,„(p), (16) 

where the subscripts c, h index the contributions by the central and the halo parts to the 
IMD given respectively by 

P lti (p) = Si(Ak = 0,K = p) (17) 

for i = c,h. Note that the IMD as expressed by eq. (|l6|) includes the possibility that the 
IMD for the bosons of the halo is different from the bosons emerging from the central part. 
Thus the relative contribution of the halo and the centre is a function of the momentum in 
this model. 

The normalization condition according to Assumption 1 reads as 

dp 

— PiAp) =1 (18) 

E 

for i = c,h. At this point the normalization conditions described for the hydrodynamic and 
Wigner emission functions can be utilized. According to Assumption 1, the hydrodynamic 
emission-function is a sum, too: 

Shyd(x,p) = (n) f c S c (x,p) + (n) (1 - f c ) S h (x,p). (19) 



It is convenient to define the 'hydrodynamic' emission function for the core and the halo, 

S h yd,c{x,v) = (n) f c S c (x,p), (20) 

S hy d,h( x 'P) = ( n ) (! ~ fc) S h (x,p). (21) 

Since these functions are normalized to the mean multiplicity of the core and the halo, 
respectively, it follows that 

(n c ) = f c (n), (22) 

(n h ) = (l-f c )(n). (23) 

Thus the core fraction f c turns out to be the the ratio of the mean multiplicity of the core 
to the mean multiplicity of the whole system. 

Integrating eqs. (p~9| - pTD over the space-time variables and using eq. (|ll|), one obtains 
that 

N 1 (p) = N 1>c (p) + N 1 , h (p), (24) 

and as a consequence of eq. (|3|) one finds that 

da da c doh 

E— = E + E . (25) 

dp dp dp 

In the forthcoming it shall be shown that the above cross-sections, describing particle pro- 
duction cross-sections from the core and the halo, can be related to measurable quantities. 
The 'ideal' shape for the BECF is expressed by 

\f c S c (Ak,K) + (l-f c )S h (Ak,K)\ 2 
C(K, Ak) = 1 + , (26) 

| f c S c (Ak = 0,K = p) + (l- f c ) S h (Ak = 0,K = p)\ 2 

which includes interference terms for boson pairs of (c, c) (c, h) and (h, h) type. Due to 
the assumption that the emission is completely chaotic, the exact value of the BECF at 
| Ak |= is always 2 in this model. 
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The measured two-particle BECF is determined for | Ak |> Q m in ~ 10 MeV/c, and any 
structure within the | Ak |< Qmin region is not resolved. However, the (c, h) and (h, h) type 
boson pairs create a narrow peak in the BECF exactly in this Ak region according to eq. (p6|), 
which cannot be resolved according to Assumption 2. From this assumption, Sh{Ak, K) ~ 
for | Ak |> Qmin or for | Ak° \> Q m in- With other words, the Fourier-transformed emission 
function of the halo for non-zero relative momenta vanishes at the given resolution Q m in- At 
zero relative momentum the same quantity gives the single particle momentum distribution 
from the halo, which is not effected by the £u>o-particle resolution. 

Including the finite resolution effect, symbolized by the horizontal bar over the correlation 



function, the measured BECF can be written as 



C(Ak, K) = l + \*(K = p; Q min ) R c (Ak, K), 

where the effective intercept parameter and the correlator of the core are defined as 

\f c S c (Ak = 0,K = p) | 2 



(27) 



R c (Ak,K) 



f c S c (Ak = 0,K = p) + (l- f c ) S h (Ak = 0,K = p)\ 2 

\S c (Ak,K)\ 2 



(28) 



(29) 



| S c {Ak = 0, K = p) | 2 

This expression indicates that the only effect of a large resonance halo is to introduce a 
momentum-dependent effective intercept parameter to the two-particle correlation function, 
while the relative momentum dependent factor shall coincide with the correlator R c (Ak, K) 
of the core. The effective intercept parameter can in turn be expressed with the IMD-P-s as 



K(K = p;Q r 



It 



A 2 (P) 



(30) 



Alternatively, the effective intercept \*{K = p\ Qmin) can be rewritten with the help of eq. 
(f22|) as 



K(K = p;Q r . 



NiAp) 

iVi(p) 



da r 



dp 



(-) 

K dpJ 



(31) 
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In the forthcoming, we shall utilize the notation A*(p) = \*(K = p; Q m in)- The effective 
intercept parameter A*(p) shall thus in general depend on the mean momentum of the 
observed boson pair, and we have < A*(p) < 1. The equation above shows that this 
effective A*(p) has a very simple interpretation as the square of the fraction of core particles 
to all particles emitted with momentum p. We also see that this equation makes it possible 
to measure the cross-section of the bosons produced in the central core or from the halo 
with the help of the combined use of the BECF parameter A*(p) and the cross-section of all 
particles. For the sake of completeness, we list these relations below: 

da c , da 

— = A(P) — ( 32 ) 
dp dp 

doh I da 

— = (1 - J Kij>)) — • (33) 
dp dp 

Alternatively, if a calculation predicts N c (p) from a model of the core, this prediction can 

be combined with the measured A*(p) to describe the IMD-N-s for the whole system as 

1 

N(p) = ^ = iV 1 , c (p). (34) 

v A *(p) 

Note that the A*(p) effective intercept parameter in our picture reflects the drop in the 
BECF from its exact value of X xct = 2 due to the contribution from the resonance halo, 
when measured with the finite two-particle momentum resolution Q m i n . Thus the mea- 
sured intercept, A*(p) does not coincide with the exact value of the BECF at zero relative 
momentum within this description, see Figure 1 for illustration. 

The effect of the halo is to introduce a momentum-dependent effective intercept param- 
eter to the BECF. The shape of the BECF for Q > Q m i n shall be solely determined by the 
freeze-out phase-space distribution in the central part. This central part or core is usually 
well accessible to hydrodynamic calculations, see refs. pl7||il~| -|T5|j for example. 



Note also that the presented results have in principle nothing to do with the similar 
relation obtained for the case of partially coherent, partially chaotic fields, where A* = ff nc , 

11 



and fi nc = (rii nc ) / (ritot) gives the fraction of the total multiplicity in the chaotic field, 
although the relation formally is similar. For the case of partially coherent fields, the Bose- 
Einstein correlation function contains an interference term in between the coherent and 



the chaotic fields which leads to a double- Gaussian or double-exponential structure [21 



emphasized recently in ref . [E21 . In our case, the measurable part of the BECF contains not 



three but only two terms, and it looks very much like the BECF would look like for the case 
the halo were missing i.e. f c = 1. 

Note also, that the focus is not on the reduction of the intercept of the BECF due to 
the halo (resonance decays) but more precisely on the momentum- dependence of this reduc- 
tion. The fact of the reduction has been known before ||^],§,[nj, however its momentum 
dependence was not utilized as far as we know. Within our formalism, it is straightforward 
to show that the core fraction of bosons can be estimated by 

, dp , 

f c = — ^Mp)A(p). (35) 

E 

Alternatively, f c can be expressed as 

(n c ) 
fc = , (36) 



(n) 
where 



E 
From eq.(J35|) it follows that the inequalities 



-dp 

n c ) = j—^X4p)N 1 (p). (37) 



miny%(p) < f c < maxy%(p) (38) 

are also satisfied. Please note that the above formulas are also valid in the special case of 
the constant function, A*(p) = const. 



Application. Present NA44 data |23[ for central S + Pb reactions at CERN SPS with 



200 AGeV bombarding energy show an approximately m t independent intercept parameter: 
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X 7T + = 0.56 ± 0.02 and 0.55 ± 0.02 at the quite different mean transverse momenta of 150 
MeV and 450 MeV Q, respectively. This suggests that in the mid-rapidity region, where 
the NA44 data were taken, the momentum distribution of pions from the resonance halo is 
similar to that from the central part, and the halo contains 1 — v'AT = 25 ± 2 % of all the 
pions. From eq. ([!(]) we can see that this observed constancy of A* in the NA44 acceptance 
regionP] leads to the following simple equalities for the IMD and the BECF: 

d 2 n d 2 n c d 2 rih 
= = , (39) 



dy dm t dy dm t dy dm 



C{Ak, K) = 1 + f 2 R c (Ak, K). (40) 

I.e. the only apparent effect of the halo is to reduce the intercept parameter of the measured 
BECF to A, = f 2 while the IMD and the (Afc, K) dependence of the BECF is determined 
by the central part exclusively. 

The resonance halo was expected to play a more active role: Pions coming from the 
resonance-halo surrounding a strongly interacting centre are predominant at lower values of 
the transverse momentum according to the SPACER calculations in ref. [|Il]] and RQMD 
calculations in ref. |T9| , similarly to the results of calculations with HYLANDER |JJ and 



the hydrodynamic calculations of the Regensburg group p4J . The resonance fraction as a 



function of the transverse momentum was explicitly shown in the publications [[0|,[| . Each 



t Note that with increasing (pt) the NA44 acceptance moves slightly in rapidity, too. For pions in 
NA44 acceptance {pt) = 150 MeV corresponds to (y) = 3.5 while {pt) = 450 MeV corresponds to 
(y) = 2.7 both values being within 0.5 units from mid-rapidity. That's why we expect the change 
in mt be the dominant effect. 

* The observed momentum-independence of the parameter X(y.mt) is further supported by the 
NA35 collaboration (which has a wide acceptance both in y and mt)- NA35 collaboration also 
finds an approximately y and mt independent \(y, mt) parameter too, as can be seen on Figs. 3 
a) and b) in T. Alber et al, Z. Phys. C66 (1995) 77. 
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of the mentioned models predicts that the halo of resonances produces predominantly low 
transverse momentum pions, and they predict that the effective intercept A*(p; Q m in) should 
increase with increasing transverse mass in contrast to the constancy of the measured values. 

The constancy of the effective parameter \*(y,mt) with respect to mt suggests a mech- 
anism of enhancing the low momentum particles in the core. Such an enhancement has a 
natural explanation in a hydrodynamic description of an expanding core as we have presented 

18| , further elaborated in a forthcoming paper |K|. The mechanism is a simple volume 



in 



effect, see also [ 15| . The calculated cross-section is proportional to the volume where the par- 



ticles with a given momentum are emitted from. This (effective and momentum-dependent) 
volume is measured by the BECF. In the NA44 experiment all three radii of this volume are 
measured and found to be equal within the errors of the measurement. Further, all the three 



radius components are found to be inversely proportional the the square root of m t |23| . This 
will give an effective volume factor, V* oc (m t ) _3//2 , enhancing the single particle cross-section 
of the core at low m t . This is to be contrasted to the case of a static fire-ball, which has a 
cross-section proportional to V, the (momentum- independent) volume of the fire-ball. This 
effect is illustrated on Figure 2. Note that the enhancement at low transverse momentum 
(as compared to the static source) is a factor of 10 larger for pions than kaons or omega 
mesons, due to the relatively small rest mass of pions. In case of pions, the concave shape 
of the momentum distribution is clearly seen on Figure 2, for heavier particles the concave 
shape is less detectable. To our best knowledge, this transverse mass dependent effective 
volume factor V* oc (m t )~ 3 / 2 , which results in a factor of 10 enhancement of low-p^ pions 
from the core, as compared to low-p* heavy resonances from the core, has not yet been taken 
into account in any previous numerical analysis of low-transverse momentum enhancement. 
Discussion. The general result for the BECF of systems with large halo coincides with 
the most frequently applied phenomenological parameterizations of the BECF in high energy 
heavy ion as well as in high energy particle reactions . Previously, this form has received 
a lot of criticism from the theoretical side, claiming that it is in disagreement with quantum 



statistics |2!| or that the A parameter is just a kind of fudge parameter, a measure of our 

14 



ignorance, which has been introduced to make theoretical predictions comparable to data. 
Now one can see that this type of parameterizations can be derived with a standard inclusion 
of quantum statistical effects - if we assume that we discuss interferometry for systems with 
large halo. Most of the studied reactions including e + + e~ annihilations various lepton- 
hadron and hadron-hadron reactions nucleon-nucleus and nucleus - nucleus collisions are 
phenomenologically well describable |2| with variations of eq. (12). This fact does not 
exclude the possibility that each boson emitting system created in the above mentioned 
high energy reactions contains a large halo. 

Let us remark that there is a standard experimental procedure to check the sensitivity 
of the fitted Bose-Einstein correlation function parameters to the variation of the bin-size. 
The data point in the first bin is frequently left out from the fitting since this data point is 
unreliable @,^3[]. These add up to a check of the stability of the fitted A*(y, m t ) parameter (s) 
for the variation of Q m in, the size of the excluded region. Clear separation of the emission 
function into a core and a halo part is possible in a given reaction only if this or a similar 
experimental test indicates the insensitivity of the fitted parameters to the exact value of 
Qmin, which may vary in a certain small relative momentum region around 10 MeV/c. 

The applicability of the halo picture to a given reaction is not necessarily the only possible 
explanation of a reduced intercept. It is known that the final state interactions may have 



an influence on the effective intercept of the two-particle correlation functions |25| and the 
effect has been shown to result in a stronger drop of the intercept value for smaller source 
i.e. for hadronic strings created in lepton-lepton, lepton-hadron or hadron-hadron collisions. 
However, more detailed calculations indicate that the higher order corrections for the final 
state interactions may very well cancel the first order effects resulting in a very small total 
final state interaction correction |[26| . In our calculations these effects have been neglected, 
so our results should in principle be compared to data which present the genuine Bose- 
Einstein correlation function. This function is not necessarily the same as the short-range 
part of the two-particle correlation function since Coulomb and Yukawa interactions as well 
as the rj'rj decay chain and other short-range correlations may mask the quantum statistical 
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correlation functions. However, as the particle density is increased (in the limit of very large 
energy or large colliding nuclei), the Bose-Einstein correlations together with the final state 
interactions dominate over other short-range correlations due to combinatoric reasons. 
It is well known that the hydrodynamic predictions fl2T| , ITfl for the rapidity and the 



transverse mass dependence of the longitudinal radius parameters of the Bose-Einstein cor- 
relation functions are found to describe the available data in heavy ion reactions at CERN 
SPS to the best available precision |^8[^3j, although the calculations were performed within 



a hydrodynamic framework, and the effects caused by resonance decays were neglected in 
the calculation. The intercept parameter was Xth = 1 in the calculation, while the intercept 
parameter in the NA35 and NA44 experiments was significantly lower than that. Thus the 
scenario discussed in this paper supports such a possibility, since it has been shown that a 
large halo of long-lived resonances results in i) a momentum- dependent effective intercept 
parameter A*(p) and ii) a special relative momentum dependence of the measurable part of 
the Bose-Einstein correlation function, which is prescribed by the emission function of the 
core. 

In summary, we have studied the case when the central boson-emitting region is sur- 
rounded by a large halo, which also emits bosons. If the size of the halo is so large that it 
cannot be resolved in Bose-Einstein correlation measurements, lot of information shall be 
concentrated in the momentum dependence of the intercept parameter of the correlation 
function. We have shown that with the help of the Bose-Einstein correlation measurement, 
the invariant momentum distribution can be measured for the two independent components 
belonging to the core and the halo, respectively. The results do not depend on any particular 
parametrization of the core nor of the halo. 

Analysis of the NA44 data for two-pion correlations indicated that the normalized in- 
variant momentum distribution of the pions from the halo of long-lived resonances within 
errors coincides with the normalized invariant momentum distribution of the pions from the 
central core. This result can be explained with a hydrodynamic model of the core evolution. 
The number of pions coming from the halo region, characterized by large length-scales of 
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Rh > 20 fm, was found to be 25 ± 2 % of the total number of pions within the NA44 
acceptance. 
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Figure captions 

Figure 1. The shape of the BECF is illustrated for a source containing a core and 
a large halo. The contribution from the halo is restricted to the shaded area, while the 
shape of the BECF outside this interval is determined completely by the contribution of 
the core. If the resolution for a given experiment is restricted to Q > 10 MeV, then an 
effective and momentum dependent intercept parameter, \*(y,m t ) will be measured, which 
can be combined with the measured momentum distribution to determine the the momentum 
distribution of the particles emitted directly from the core. 

Figure 2. The momentum distributions d 2 N/(mf dydmt) are shown for pions, kaons and 
omega mesons in arbitrary units as a function of m t — m, when the effective m« dependent 
volume factor, V* oc l/m t is taken into account. To illustrate the effect of the volume 
factor V* an effective temperature T e ff = 200 MeV was used for each particle, and we 

utilized m n = 140 MeV, mx = 494 MeV and m^ = 782 MeV. The plotted functions are 

1 
all proportional to exp(— (m t — m)/T*), being normalized to have the same value at 

(3/2) 

m t 
m t — m = 2.0 GeV. Note that the enhancement at low transverse momentum is a factor of 

10 larger for pions (solid line) emitted from the core than the enhancement for the heavier 

omega mesons (dashed line). The kaons (short-dashed line) and omega mesons on the 

other hand have similar momentum distributions. Dotted line stands for the momentum 

distribution of a static source, which is proportional to exp(— m t /T*). 
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